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Abstract 

We present a model in which metastable supercooled phase and stable equi- 
librium phase of vortex matter coexist in different regions of a sample. Minor 
hysteresis loops are calculated with the simple assumption of the two phases of 
vortex matter having field-independent critical current densities. We use our 
earlier published ideas that the free energy barrier separating the metastable 
and stable phases reduces as the magnetic induction moves farther from the 
first order phase transition line, and that metastable to stable transformations 
occur in local regions of the sample when the local energy dissipation exceeds 
a critical value. Previously reported anomalous features in minor hysteresis 
loops are reproduced, and calculated field profiles are presented. 
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I. INTRODUCTION 



Supercooled or metastable states have been reported across first order phase transitions 
in vortex matter |l -10|. It has also been established that isothermal field excursions cause 



the metastable phase to be converted to the stable phase []5|,|7|,|8|.|T0[| . We have proposed (TT 



that the isothermal field variations provide a fluctuation energy that causes the metastable 
supercooled phase to cross the free-energy barrier and transform to the stable equilibrium 
phase. 

The experimental techniques used to study the magnetic signatures of such transforma- 
tions in vortex matter are (i) bulk dc measurements using a SQUID or a vibrating sample 
magnetometer which yields the magnetisation M of the entire sample PHlQl; (ii) bulk ac 
measurements of susceptibility which probe a region near the surface of the sample 0; and 
(iii) local measurements of magnetic induction using magneto-optic or microhall probes that 
allow mapping the spatial profile B(x) |1]|||4|]. The last of these three techniques has been 
used recently to show that different (metastable and equilibrium) phases exist simultane- 
ously in different regions of the sample . Similar inference has also been drawn by studies 
using the dc magnetisation |5|,[| and ac susceptibility || techniques. Experiments have thus 
shown that metastable to stable transformations occur over local regions, and prompted by 
these developments we have given a formalism to calculate spatially-resolved energy dissi- 
pation under an isothermal field variation [Q. In this paper we shall use this formalism to 
calculate how the metastable to stable transformation progresses inwards from the surface 
of the sample, under experimentally relevant field excursions. We shall also calculate the 
sample magnetisation M and the spatial field profile B(x) as the experimental observables. 

II. MODELLING THE PEAK EFFECT 

The response of a hard superconductor to external magnetic fields is understood in 
terms of Bean's critical state model (CSM) ||13|| . Bean had assumed that the critical current 
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density Jc is independent of field. While detailed agreement with experiments has required 
introduction of various functional forms of Jc(B) (see e.g. ref. [14]), much of the essential 
physics is captured even by assuming a field-independent Jc- 

Experiments have recently been addressing the region below and near the onset of a 
peak in Jc(B) at B=Bx, where a first order phase transition is seen in some superconductors 
|lHT0|. The occurrence of this "peak-effect" has been known in various superconductors for 
a very long time, but attempts to have a CSM describing this Jc(B) have been made only 
recently While our detailed analytical model |L5j could be used for the subsequent 



calculations, our focus here is to understand whether qualitatively new and anomalous signa- 
tures in recent experiments p|- |T0|l can be arising from metastable to stable transformations 
in vortex matter. In this paper we shall use a simple Bean-like assumption for the two 
phases of vortex matter viz. 

J C (B) = Ji for B<B 1 

= J 2 for B > B 1 (1) 

We stress that, at a fixed temperature, our model has only two constant parameters viz. 
(J2/J1) and Bi. We have assumed above that phase 1 (characterised by Ji) is the stable 
phase for B(x)< B%, and phase 2 (characterised by J2) is the stable phase for B(x)> B\. We 
recognise that Jc is not a thermodynamic quantity, but is a physical property that changes 
disco ntinuously across the phase transition. Phase 2 can exist for B(x) < B\ as a supercooled 
metastable phase, and we shall address this possibility in the next section. In this section 
we shall assume, however, that the free energy barrier surrounding the metastable phase 
drops very sharply as B(x)=Bi is crossed, and supercooling or superheating does not occur. 

To obtain magnetisation-vs-field (or M-H) curves, we consider the sample to be in the 
form of an infinite slab in parallel field, as this geometry has the simplest algebra amongst 
the zero demagnetization factor cases of infinite cylinders in parallel field. We shall also 
continue with Bean's simplifying assumption of Hci=0 followed usually in the CSM Ji"3]-|r5 



We follow standard procedures ]14 , |15| to solve the CSM, and show in figure 1 the 



envelope M-H curves obtained with equation (1), with the parameters Bi = 1500 mTesla, 
J 1 R=4mTesla, and J2R=10 mTesla. Here the slab has surfaces at x=±R, is infinite along 
the y and z directions, and the magnetic field is applied along the z-axis. (We shall consider 
only positive values of x in this paper; there is a symmetry about x=0.) Note that the first 
order transition shows different widths in M vs H when measured along the field- increasing 
and along the field-decreasing directions. This is because the shielding current density at x 
is dictated []T3| |3~^1 by the local magnetic induction B(x) through equation (1), and B(x) is 
different from the applied field as well as different in the field-increasing and field- decreasing 
cases. In figure 2 we plot B(x) for some values of applied field H corresponding to the field- 
increasing and field-decreasing cases. We note that phase 1 and 2 exist simultaneously in 
two different regions of the sample. We emphasize that there is no metastability because the 
stable phase 1 exists wherever B(x)< Bi and the stable phase 2 exists wherever B(x)> B\. 

The calculation above is for some fixed temperature Ti, and we note that the phase 
transition field Bx falls as the temperature Ti rises 



III. SUPERCOOLING AND METASTABLE-TO-STABLE TRANSFORMATIONS 

We now consider that we have applied a field Hi which is smaller than Bi(Ti). But we 
apply this field at a much higher temperature T 2 such that Hx is much larger than B^T^. 
So, B(x) throughout the sample is larger than Bi at that temperature, and the entire sample 
is in phase 2. We assume further that B(x) is constant at T 2 . This happens if the critical 
current density J 2 in phase 2 vanishes at T 2 . One can, however, also achieve a constant 
B(x) by applying an external field H!+hcos(wt), with (H 1 -B 1 (T 2 )) > h > J 2 (T 2 )R, and 
then slowly reducing the amplitude h to zero p^Jl8[] . 

We now lower the sample temperature (i.e. field-cool) to Ti such that the sample is 
supercooled and is metastable in phase 2. As discussed in references [11,12,17], there is a 
free energy barrier fe(T) that keeps phase 2 metastable, where fs(T) is determined uniquely 
by B(x) and T. The vortex matter in the neighbourhood of x will transform to phase 1 when 
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the fluctuation energy Pd(x) created by an isothermal field variation is larger than [fe(T) - 
kT]. 

The field profile B(x) in the field-cooled sample is constant at Hi (see figure 3), and we 
now start lowering the applied field H with the temperature fixed at Ti. Since the sample is 
in the supercooled phase 2, the shielding currents set up initially will have a magnitude J2. 
The variations in field will cause a fluctuation energy P<f(x) given by equations (4) and (5) 
of reference [12], and the vortex matter in the neighbourhood of x=xo will transform to the 
stable phase at Pd(x )= £b(T) - kT =Pq- This transformation is triggered from the surface 



1 1^| and the shielding current magnitude will drop to Ji for x > Xq. The point xo moves 
from (x/R)=l to (x/R)=0 as the applied field is lowered, and B(x) are shown in figure 3 
for representative values of the applied field. We have used Hi=1480 mTesla, and Po=18 
(mTesla) 2 . The large (small) slopes of B(x) correspond to large (small) magnitudes of the 
shielding current density, and thus to vortex matter being in phase 2 (phase 1). From these 
B(x) one can readily calculate p3|-4T5|1 the sample magnetisation as the field H is lowered. 
In figure 4 we show the minor hysteresis loop (MHL) obtained as the applied field is lowered 
after field-cooling. We show also the field-decreasing envelope curve from figure 1. Note 
that the MHL first shoots out above the envelope curve, and then slowly merges from above. 
This nature is in qualitative agreement with published data |7|j9|, ^0| , pj| , p0[| . If we had used 
the detailed Jc*(B) of reference [15] to model the peak-effect, instead of the simple model of 
equation (1), the peak of the MHL would be less sharp and the merger with the envelope 
curve would be slower. The simple model used to obtain figure 4 brings out the qualitative 
behaviour observed and captures the essential underlying origin of anomalous MHLs as being 
due to phase 2 being supercooled and the transformation from the metastable phase 2 to 
the stable phase 1 occurring progressively deeper into the sample. 

We have assumed that P (and thus fs(Ti)) is constant at 18 (mTesla) 2 . fs(Ti) is 
actually dictated by B(xo), and falls monotonically as B(xo) moves farther from the phase 
transition line Bi(T). As is seen in figure 3, B(x ) varies only by less than a few mTesla as 
the MHL merges with the envelope curve. The assumption of a constant Pq over an MHL is 



thus justified. If, however, we field-cool to the same temperature Ti at a lower field H 2 , then 
fe(T) will be lower P77[] . This implies that Po(Hfc=H 2 ) will be smaller than P (H F c=Ki). 
We show, in figure 5(a), the MHL for the case when the sample was field-cooled to H2=1400 
mTesla where Po is taken to be 4.5 (mTesla) 2 . In figure 5(b) we have taken 13^(7=153=1300 
mTesla, where fe must be still lower and is taken as Po=2 (mTesla) 2 . The MHLs again 
shoot out of the envelope curve, but to peak values progressively smaller than in figure 4. 
The merger of the MHLs with the envelope curve also occurs over a progressively narrower 
range of field reduction than in figure 4. This qualitative change in the nature of the MHLs 
with reduction of B.pc is also consistent with published data |7|fl p!0| , p^ .^0|l . 



IV. CONCLUSION 

We have used the ideas developed in references [11,12,17] to calculate the isothermal 
field-cooled MHLs, and the spatial field profiles B(x). The model calculation was done, in 
the spirit of Bean's original work with field-independent critical current densities. The 
only parameters were (J2/J1), and the onset field Bi at which the peak effect starts in the 
field increasing case. We used the fact [17] that fg becomes smaller as B(x) falls below Bi, 



and that metastable to stable transformations occur in local regions of the sample [12 



The formalism of reference [12] can similarly be used to calculate MHLs after different 
thermomagnetic histories. We assert here that the simple model of equation (1) reproduces 
qualitative features of various observations [p|- [T0| , p!9|j20|] of anomalous MHLs. As was stated 
in the Introduction, more detailed tests of the extent of phase coexistence are possible and 
calculated B(x) can be compared with field profiles measured with local probes. Our model 
also predicts the spatial region over which the two phases coexist, and the evolution of these 
regions under isothermal field variation. 

We thank Drs. S.B. Roy and Sujeet Chaudhary for many helpful discussions. 
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FIGURES 

FIG. 1. Field-increasing and field-decreasing envelope M-H curves are shown, following the 
model of equation (1), with Bi=1500 mTesla. The crosses indicate applied field values at which 
B(x) profiles are shown in figure 2. 

FIG. 2. B(x) profiles are shown for (a) field-increasing case at applied fields of 1495, 1505, 
and 1515 mTesla; and (b) field-decreasing case at applied fields of 1501, 1498, and 1495 mTesla. 
The arrows indicate the x at which B(x) = 1500 mTesla. 

FIG. 3. The sample is field-cooled in 1480 mTesla, when B(x) is constant as shown by the 
dashed line, and vortex matter is in the metastable phase 2. B(x) are shown as the applied field is 
lowered isothermally to 1475, 1473, and 1471 mTesla. In the last two fields the vortex matter has 
transformed to the stable phase 1 at x>xo, where xo is indicated by an arrow. 

FIG. 4. MHL obtained after field-cooling at 1480 mTesla is shown by the solid line. It 
overshoots the envelope curve, and merges with it slowly from above. The circle indicates the 
starting point of the MHL, with M=0 corresponding to the constant B(x). 

FIG. 5. Same as figure 4 except that field-cooling was done at (a) 1400 mTesla; and (b) 1300 
mTesla. 
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